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Capillary origami: spontaneous wrapping of a droplet with an elastic sheet 
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Abstract The interaction between elasticity and capillarity is used to produce three dimensional 
structures, through the wrapping of a liquid droplet by a planar sheet. The final encapsulated 
3D shape is controlled by tayloring the initial geometry of the flat membrane. A 2D model shows 
the evolution of open sheets to closed structures and predicts a critical length scale below which 
encapsulation cannot occur, which is verified experimentally. This elastocapillary length is found to 
depend on the thickness as /i 3 ^ 2 , a scaling favorable to miniaturization which suggests a new way of 
mass production of 3D micro- or nano-scale objects. 

PACS numbers: 46.32,+x, 68.08.-p, 81.16.Dn, 85.85.+j 
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Origami and haute couture design both rely on folding 
and assembling planar material to create elegant three- 
dimensional (3D) shapes whose variety and complexity is 
governed by the number, order and orientation of folds. 
Folding is also a way to reduce the size of deployable 
structures in space industry (satellite pannels or sun 
sails [? ]) as well as in nature (plant leaves folded in 
a bud [? ]). In micro and nano-fabrication, the fold- 
ing of planar structures is a promising approach to build 
3D objects, since most microfabrication technologies pro- 
duce planar layers through surface etching [? ] . 

While centimeter-sized objects may be folded using 
well disposed magnets [? ], capillarity is a particularly 
relevant mechanism for micro-fabrication [? ], as sur- 
face forces dominate over bulk forces at small scales. In- 
deed, the capillary attraction and stiction of wet slender 
structures may induce disastrous damage in micro and 
nano-scale fabrication [? ? ? ? ], in addition to being 
involved in lung airway closure (neonatal respiratory dis- 
tress syndrome) [? ] . Capillary interactions have already 
been proposed as a way of assembling and orienting rigid 
objects in 2D at the surface of water [? ], or in 3D by cap- 
illary bridges [? ], while further studies have focused on 
the capillary induced deflection of elastic rods [???]. 
Recently, the multi-folding of a flexible ribbon squeezed 
in a meniscus has also been described through an elegant 
experiment [? ]. Here, we address the effect of capillary 
forces on a flat elastic membrane and we show how the 
spontaneous folding of a flexible sheet around a liquid 
droplet leads to a predetermined 3D shape. 

Our experiments were conducted using polydimethyl- 
siloxane (PDMS) membranes. The PDMS (Dow corning 
Sylgard 184, 10:1 polymer /curing agent mix) was spin 
coated at 24° C on a glass microscope slide, at rotation 
rates of 1000-2000 rpm. Once the PDMS was cured, this 
resulted in sheets with thickness in the range 80-40 /im. 

During a typical experiment, a geometric shape is man- 
ually cut out from the PDMS layer and placed on a super- 
hydrophobic surface. A drop of water, of volume 1— 80 /iL 
depending on membrane size, is deposited on the PDMS, 




FIG. 1: Wrapping of a drop of water with (a) square and 
(b) triangular PDMS sheets. The total elapsed time for each 
experiment is about half an hour. 



while making sure that the water reaches all the corners 
of the flat sheet. The water is then allowed to evapo- 
rate at room temperature and a time lapse sequence of 
images is taken of the membrane/drop pair. As the wa- 
ter volume decreases from its initial value to complete 
evaporation, the surface tension of the liquid pulls the 
sheet around smaller volumes, thus increasingly curving 
it (Fig. [I]). Sufficiently thin sheets eventually encapsulate 
the liquid with a shape that depends on the initial cut of 
the membrane. 
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In the simplest drop is deposited on a square 

PDMS sheet, as shown in Fig. This initially leads 
to the bending of the four corners towards the center 
(mode-4). As the drop volume decreases further, this 
bending mode becomes unstable and the shape rapidly 
switches to a mode-2 state; in this state, the corners are 
attracted to each other two by two and the flat central 
region disappears, giving rise to a quasi-cylindrical shape 
(middle image in FigQJi). As the drop volume decreases 
further, the approaching edges eventually touch and a 
water-filled tube, a few hundred microns in diameter, is 
formed. Finally, the water evaporates completely and 
the tube elarges into a loop whose ends merge into a flat 
contacting region, similar to the shape shown in Fig. [2^.. 

A different situation occurs when a drop is de- 
posited on a triangular sheet, see Fig. [TJd and Movie SI 
(www.pmmh.espci.fr/~benoit/wrapping.avi). The cor- 
ners fold towards the center as in the previous case, but 
this leads to a stable mode-3 folding. The three cor- 
ners eventually join at the center and the sheet seals into 
a tetrahedral pyramid. At this stage, the curvature of 
the surface is concentrated on the edges of the pyramid, 
while the sides are almost flat (fourth image in Fig. [TJd). 
Further evaporation, which now occurs at a slower rate, 
finally leads to the buckling of the pyramid walls due to 
the negative internal pressure (final image). 

For both of the above geometries, the scenario is dif- 
ferent for membranes of smaller size or larger thickness. 
As the stiffness increases, the forces generated by sur- 
face tension become insufficient to achieve complete clo- 
sure: after going through similar initial stages as for thin 
PDMS layers, a maximum bending of the sheet is pro- 
duced after which it reopens again, never having pro- 
duced an encapsulated drop. The final state is there- 
fore a planar sheet covered by a thin film of liquid which 
quickly evaporates. 

The membrane deformations reduce the liquid-air area 
A and thus the surface energy jA (7 being the surface 
tension), at the cost of increasing the elastic energy. For 
a thin plate, the isometric bending energy density is de- 
fined locally as as Bn 2 /2, where n is the curvature and 
B = Eh 3 /12(1 - v 2 ) is the bending stiffness (E is the 
Young's modulus, v is Poisson's ratio and h is the thick- 
ness) [? ]. For a system of typical dimension L, folding 
will occur if the surface energy (jL 2 ) is large compared 
to the total bending energy (B). This is equivalent to 
comparing L to 
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that we refer to as the elastocapillary length [? ]. Physi- 
cally, L ec gives the typical radius of curvature generated 
by capillary forces for a given bending rigidity. 

The elastocapillary length can be measured directly 
through a calibration experiment rather than by estimat- 
ing each parameter in Eq. [1] which would produce large 
errors. The experiment consists in depositing a drop of 
wetting liquid onto a very long strip of PDMS, which 



folds over into a loop with a self-contacting tail. We 
consider the late stage, when only a small meniscus still 
holds the loop together (Fig. [5^). 
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FIG. 2: (a) Experimental image of a loop of PDMS held by 
a meniscus of ethanol, and comparison with the rescaled the- 
oretical shape (white curve), (b) Virtual displacement dx of 
the meniscus position. 



This equilibrium state is 2D, so the following analysis 
is understood to be per unit depth. The cross-section 
follows planar rod elasticity coupled with surface tension. 
We note 6{s) the angle made by the tangent to the rod t 
with respect to the horizontal, at curvilinear coordinate 
s. We also note R the constant vectorial tension of the 
beam, so that Euler's elastica equations for equilibrium 
may be expressed as [? ] 



B— e z + t xR = 0, 



(2) 



where e z is the unit vector perpendicular to the plane of 
the section. The global shape of the loop is determined 
by elasticity, given the curvature at the contact point. 
Since this curvature is set by an equilibrium between sur- 
face tension and elasticity, it must scale as 1/Lec- The 
prefactor is found by considering a virtual displacement 
dx of the contact point for a loop of curvilinear length 
2£ (see Fig. [3d). The elastic energy, E ei = f_ £ Bn 2 /2ds, 
varies due to the change in the integral bounds, the inte- 
grand being stationary near elastic equilibrium. It may 
be written as dE e \ = Bn 2 dx, where n c is the curvature 
at the contact point. On the other hand, the variation 
of surface energy due to replacing a solid-gas interface 
(surface energy jsg) with a solid-liquid interface (jsl) 
is 2(7sx — jsG)dx = 2j cos adx, where a is the contact 
angle of the liquid on the surface (here we use ethanol 
which has a zero contact angle a on the PDMS). Bal- 
ancing the two energies yields the value of equilibrium 
curvature at the contact point: 
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The boundary value problem of Eqs. ([2]) and (j3]) was 
solved numerically by enforcing the contact condition at 
the base of the loop, J^_ £ sm(0)ds = 0, and Lec was mea- 
sured by scaling the numerical solution to fit the experi- 
mental shape. An excellent agreement is obtained (white 



3 



line superimposed on experimental picture in fig[2^). In 
practice the width of the loop (theoretical value 0.89Lec) 
is sufficient to determine Lec- 
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FIG. 3: Critical length for folding vs. elasto- capillary length: 
black squares: square sheets, black triangles: triangular 
sheets (leading to 3D pyramids), lines: linear regressions. 



The value of L ec was measured for different membrane 
thicknesses and for each thickness, the minimum length 
for encapsulation L crit was determined for square and 
triangular membrane shapes. We found that L crit varied 
linearly with Lec for both shapes, as one might expect 
from dimensional arguments. The slopes were 7.0 for the 
squares and 11.9 for the triangles (Fig [3]). 

An understanding of the evolution of the folding is ob- 
tained from a 2D theoretical description of the elasto- 
capillary balance: let us consider an elastic inextensi- 
ble rod of length L bent by pressure and surface tension 
forces at its ends (Fig. [4] inset). Equation ([2]) still holds 
but the vectorial tension R is not constant anymore due 
to the pressure forces but follows 



ds 



pn = 0, 



(4) 



where n(s) is the unit vector normal to the rod. By 
neglecting gravity, the pressure p in the drop is uniform 
and the 2D shape of the liquid-air interface is a circular 
arc of radius r = j/p. 

Equations ([2]) and (|4]) , with appropriate boundary con- 
ditions, are solved numerically and the distance 5 be- 
tween the rod tips is calculated as a function of drop 
section area S for different rod lengths L, as shown in 
Fig.[H Two trivial flat states (S/L = 1) exist in all cases 
for an infinitely large drop (A and A') or when no liquid is 
present (B or B'). However, the evolution between these 
states is qualitatively different depending on L/Lec- 

For L/Lec slightly below the critical value (Fig. EK), 
these two states are connected by a continuous family 
of solutions: the drying of a large drop leads to tran- 
sient limited bending but eventually to reopening. Closed 



states (C) do exist in a restricted domain of the phase 
space, as the stable segment of a solution loop, while the 
unstable part of this loop corresponds to a drop which has 
depinned from the edges (D). This bistability is indeed 
observed in the experiments: while small sheets never 
close spontaneously, they can be forced into a closed state 
if the liquid is pulled out with a syringe. 

Above the critical ratio L/Lec, the drying of a large 
drop (A') leads continuously to complete wrapping (C), 
as observed experimentally in Fig. [H Another branch 
links the open state (B') to unstable solutions (D'). How- 
ever, this branch is difficult to observe experimentally 
because it requires the spreading of a thin film of water 
on the flat membrane. Finally, figure [4] shows that the 
initial conditions play a major role in selecting the final 
observed state because of the presence of bistability. 
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FIG. 4: Distance between the rod tips versus drop volume 
for different values of the rod length L: (a) L = 3ALec re- 
opening of the rod as S decreases; (b) L = ALec complete 
wrapping of the rod around the drop. The inset shows a 
sketch of a 2D drop deposited on a flexible sheet. 



In the numerical simulations, the crossing from dia- 
grams (a) to (b) was found at L cr u = 3.54 Lec- This is 
below the value obtained for the square membranes, even 
though the final state of a folded square is 2D (Fig. []Ji). 
This discrepancy is likely due to gravity and 3D effects, 
neither of which is accounted for in the model. Indeed, 
the drops considered in Fig. [3] are slightly flattened by 
gravity, since they are larger than the capillary length 
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(L c = j/pg, where p is the density and g the acceleration 
of gravity). More importantly, the early stages are always 
3D as seen in the folding of the corners (e.g. mode-4 fold- 
ing) and thus deviate the evolution from the prediction 
of the 2D model. A 3D theoretical model should include 
the full description of thin plate elasticity. Only such a 
model would account for the localization of the curvature 
around the edges of the pyramid for example, which are 
reminiscent of the crumpling of thin sheets. 

More generally, 3D wrapping of a droplet leads to a 
geometrical incompatibility described by Gauss's theo- 
rema egregium. A classic example is that it is impossible 
to draw a planar map of the earth which conserves dis- 
tances: the wrapping of a spherical object with a planar 
sheet must involve stretching as well as bending. For 
a thin membrane, the stretching energy scales as the 
thickness h and is therefore much more expensive than 
the bending energy which scales as h 3 . The unavoid- 
able stretching is therefore localized in singular regions, 
leading to crumpling singuliarities [? ? ]. We expect 
that for L ^> Lec, the number of singularities will in- 
crease since surface tension becomes dominant over elas- 
tic forces, leading to quasi-spherical encapsulated shapes. 

Precise engineering of the final closed state may be ob- 
tained by tayloring the initial sheet geometry. In this 
way, it is possible to approximate a smooth sphere by 
starting with a flower shape, as shown in Fig. [5^. Fig- 
ure [5)3 shows how to form a cube by starting from a 
cross shape. As in the case of the pyramid, evaporation 



after encapsulation takes place at a reduced rate and may 
buckle the sphere or the cube due to the negative pressure 
thus generated. However, one may also freeze (or poly- 
merize) the encapsulted liquid at any point in this pro- 
cess, therefore fixing the desired shape in place. Finally, 
small perturbations of the initial shape may also be used 
to yield different final states. In particular, rounding-off 
two opposite corners of a square sheet leads to the mem- 
brane closing along its diagonal (Fig. [St) rather than par- 
allel to the sides. These experiments suggest that a wide 
variety of final shapes may be achieved through careful 
tuning of the initial 2D shape. 

In summary, three dimensional sub-millimetric objects 
are produced from 2D membranes through the interac- 
tion of elasticity and capillarity. The limitation on the 
minimal size for folding is determined by the elastocap- 
illary length (Eq. [1]) which scales as the sheet thickness 
^3/2 This scaling is favorable to miniaturization since 
thinner membranes lead to much smaller critical lengths. 
This opens the way to mass production of micron-scale 
3D structures based on standard microfabrication meth- 
ods since a wide variety of objects can be fabricated by 
tayloring the initial planar sheet. 
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FIG. 5: Tuning the initial shape to obtain (a) a spherical 
encapsulation, (b) a cubic encapsulation or (c) a triangular 
mode- 2 fold. 



